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ABSTRACT 

A   simple    theoretical    relation    between   a      one- dimensional    sur- 
face  wave  and   the   amplitude    fluctuations    of  an    acoustic    standing 
wave    in   a    cavity   with    idealized  boundaries    is   presented.       An    ex- 
periment  was    designed   to    investigate    the  applicability    of    the 
theory    to  a   more   complicated  configuration.      A   237    cm  x    117   cm 
x   9    cm  water-filled  cavity   with  pressure-release   walls   was    excited 
in   a    normal    mode   by   a    source  placed    in    one   corner.      The   resulting 
acoustic   signal    was   monitored  at    various    locations    throughout    the 
cavity   with   a    small  probe   receiver.       Wind-driven    surface   waves 
were   generated   by   a   multi-fan   blower    placed  at    one    end    of    the 
cavity.      The   relationship    between    the   resulting    fluctuation    of   ths 
acoustic    signal   and   the    surface   fluctuation   above   the   receiver    was 
investigated.       The   results    are   consistent   with    the   simple    theory 
at    low   surface-wave   frequencies.      Anomalous    behavior    was    observed 
for    the  higher    frequencies    of   the   surface  wave   spectrum.      Possible 
mechanisms    for    this    inconsistency   are  postulated  and   discussed. 
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I.   INTRODUCTION 

In  A.B.  Wood's  model  study  of  shallow  water  propagation  LlJ,  a 
sound  amplitude  fluctuation,  synchronized  with  the  passage  of  sur- 
face waves,  is  reported.   Wood  suggested  the  cause  was  an  appre- 
ciable change  in  water  depth  in  the  vicinity  of  the  source  and/or 
specular  scatter  from  the  surface.   The  reflection  and  scattering 
of  sound  from  uneven  surfaces"  has  been  given  considerable  attention 
in  the  past  twenty  years  [2] .   Most  treatments  in  the  open  litera- 
ture have  involved  ray  theoretic  approaches,  treating  the  reflect- 
ion and  scattering  of  plane  waves  from  irregular  surfaces,  and 
have  assumed  no  interaction  with  a  second  boundary.   When  a  second 
boundary  is  present,  the  possibility  of  ducted  propagation  arises. 
A  practical  mathematical  tool  for  investigation  of  surface  irregu- 
larities affecting  this  type  of  propagation  would  be  normal  mode 
theory.   Normal  mode  approaches  L3j4j  have  been  rare. 

Several  experimental  investigations  involving  normal  mode  pro- 
pagation have  been  conducted.   J. A.  Scrimger  [ 5J  reported  corre- 
lation between  surface-wave  heights  and  sound-amplitude  fluctuat- 
ions both  in  magnitude  and  in  frequency  content.   O.S.  Tonakanov  [6] 
found  a  linear  relationship  between  surface-wave  heights  and  per- 
centage fluctuation  in  sound  amplitude  and  phase. 

The  present  work  describes  an  experimental  method  for  the  in- 
vestigation of  the  effect  of  surface  variation  on  sound-fluctuat- 
ion in  guided  mode  propagation.   Sound-amplitude  fluctuations, 


caused  by  a   wind-driven    surface,    were   measured  and    compared    to 
the    surface    irregularity    in    the   vicinity    of    the   receiver.      The 
surface   was    treated   as    a    Fourier    superposition    of    various    sinu- 
soidal  waves,    and    the   experimental    results   at   each    frequency   were 
compared   to  a    theoretical   prediction    for    sound-amplitude   fluctu- 
ations   based   on   a    monof r equency   sinusoidal   wave   surface. 
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II.   THEORY 

A  normal  mode  solution  to  the  wave  equation  for  a  rectangular 

enclosure  with  smooth  boundaries  (Figure  1)  can  be  expressed  as 

$   =  cos  k  x  cos  k  y  sin  k  z  cos  COt  (2.1) 

x         £  j 

where   CO   =  circular  frequency  of  excitation 

th 

k.    =  propagation   constant    for    the    i —  direction. 

If   the   boundaries   perpendicular    to    the   z    axis    are  pressure- 
release   and  are    located   at    z    =   0   and   z    =  I    ,    then    k^l      =   n_    is 

z  3  z  3 

required.   In  all  that  follows,  it  is  assumed  that  the  system  is 

excited  at  frequencies  allowing  only  n   =  1  so  that  the  eigenvalue 

3 

k„  is  as  small  as  possible. 
3 


Figure  1.  Theoretical  Cavity  for  Guided  Mode 
Propagation 

If  the  upper  boundary  is  disturbed,  the  boundaries  in  the  z 
direction  then  occur  at 

z  =  0 


z  =  I       +   GF(x,y,t) 


(2.2) 
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where  ep(x,y,t)  describes  the  change  in  depth  of  the  cavity  re- 
sulting from  the  surface  fluctuation. 

If  GF(x,y,t)  is  small,  then  a  solution  of  the  form 


=  V  C02 


(2.3) 


is  sought  xvhere  0..  is  the  solution  for  smooth  boundaries  shown 
above  and  G0   expresses  a  perturbation  due  to  the  altered  boundary 
at  z  =  H    . 

Z 

When   e    is    sufficiently    small,    the   following    boundary   conditions 


on   02   result   [7] 


0. 


-&0, 


z=i 


6z 


:=£ 


F(x,y,t) 


(2.4) 


0, 


=    0 


z=0 


and  0   must  satisfy  the  wave  equation 


D  tfo  =  o. 


If  the  surface  fluctuation  is  a   monofrequency  wave  propa- 
gating across  the  surface  z  =  £      in  the  x  direction,  F(x,y,t)  may 
be  written 


F(x,y,t)  =  I      cos(yx-nt) 


(2.5) 


where  y    =   surface   wave  propagation    constant 

Cl  =   circular    frequency    of   the    surface  wave. 
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The   solution    for    0      under    these   conditions    is 


02    =  t   k3ZzCOS    k2V 


s  in   (X    z 

—. ±-r-  cos[(k_-ty)x-((*>f^)t]  + 

sxn  \i,.i 


1    z 


sxn   \i    z 

sin   |i    £ 

2     2 


cost  ( k  ~y ) x+  (ohH)  t  ]  + 


sin   |i    z 

— : V~  cos[(k1+y)x+(W-fi)t]  + 

sin   |i    X  v     1       '       v 

3  z 


sin  H , z 

— : ^r-   cos[(k-y)x-(CO-ft)t] 

sxn   |x,£  v    1       ' 

4    z 


(2.6) 


where 


(^)'-(V^)2-k22 


(2.7) 


CO+Q 
*     c    ' 


(kry)2  -  k22 


*     c    ' 


2  2 

(k1+y)    -  k2 


<*£>   -  i\-y)2  -  *22 


If  €1  «   CO,    then    the   approximation 


2-         2     •    ,  2    ,,  v2,2 

^1      =  ^3      =        "(k1+>/)    "k2 

2    '  2    ■    .  2        '      „.2    ,     2 

^2     =  ^     =  k  -(kry)    -k2 


(2.8) 


can  be  made  and  the  complete  solution  simplifies  to 

$  =  cos  cot  cos  k  v[cos  k  x  sin  k0z  +  -  e k  i    Acos  (Qt+0 ) ]       (2-9) 
2        1        o  ^         o    z 


where 


■N 


\i    z    \         /sin  p<2z  \  /sin   M^z  \  /sin   M^z 

7l~T~  )    +  I  sin   [i^l     J    +   2  (sin   li^JUin   V  ^ .    '  °C         ^ 


sxn 


2    z 


2    z- 


1     Z) 


2    z 
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and 


sin   p    z   \  /sin   p-z 

1  si"(Vr>*  +(Sinu  i  J  si"(krr) 


sin   p,£ 

1     Z 


2    z 


Tan;/)    =    - 


sin   p    z 


sin   p.    2 


___     cos(k1+y)x  +     -I__     cost^-yjx 


1  2 


2    z 


Notice   that    it    is    the   second   term    in    the    brackets    of   Eq.    2.9    that 
describes    the   fluctuation    in    the   acoustic   amplitude  and   that    this 
term    depends    only    on   x   and   z.      For    the   subsequent   analysis    Eq.    2.9 
will    be  written    in    the    form 


=   cos   LOt    cos    k   y   [cos    k  x   sin    k   2    +    a   cos(f2t  +0 )  ] 


(2.10) 


where        k,^      =  n 
3  z 


c/e    =(tt/2)a. 


When    k  x    =  TT    (an   antinode    of    the   standing  wave    in   x)  , 
c/e    becomes 


c/e    = 


s  ]  n  p    z 

I  sin   p,£ 
1    z 


+    2 


fsin   p.    z 


sin   P9z 


sin   li  J      I      sin   ^J 
1    z  /   \  2    z 


sin   P2z 

sin   \h    l 
2    z 


(2.11) 


At   z    =   1/2,    this    expression    further    reduces    to 


a/e    =  tt/4 


cos    \l.Ji 
1    z 


cos    po£ 
2    z 


and 


$    =   cos   COt    cos    k   y    [l  +  c   cos(S7t  +  0)] 


(2.12) 


(2.13) 
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Predicted  Values  of  a/c 
vs 
Surface  Wave  Frequency  (^/2tt) 
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On    the   other    hand,    at    a    node  where   k  x   =  rr/2   and  for 
x   =    I    /2,    Eq.    2.12   becomes 


o/e    =  tt/4 


cos   u   £ 
1    z 


COS     \l     I 

2   z 


(2.14) 


and        $    =   a   cos    OJt    cos    k   y    cos(Qt  +  0)  .  (2.15) 

Plots    of   a/e    as    functions    of   surface   wave   frequency   for    the 
two   cases    are  presented    in    Figures    2   and    3-      These    functions    dis- 
play   strong    dependence    on   both   acoustic   and   surface-wave  propa- 
gation   constants. 

Notice    that  when  y    -   k    ,    there    is   a    "resonance"   effect.      The 
perturbation   analysis    may   fail    in    this    region,    because   C0      may  no 
longer    be   small    compared   to  0    .       If   bulk  absorption    of    the   medium 
is    included,    a/e    no    longer    becomes    infinite,    but    is    still    extremely 
large.      For    y  <   k    ,    the   surface   fluctuation    tends    to   behave  as   a 
rather    uniform   change    in    depth    of    the   cavity,    and   the   amplitude 
fluctuation    behaves   as    i'f   the   standing  wave    in    the   z    direction 
were   squeezed  and  then    released.      When  y  >    k    ,    the   amplitude 
fluctuation   appears    to   react    to   the   average   fluctuation    over   many 
surface   wavelengths,    and   approaches    the  steady-state   solution   0    . 

The  relation    between    frequency   and  propagation    constant    is 

CO/k   =   c  (2.16) 

for    the  acoustic   wave   and 

tyy  -  v  (2.17) 

for  the  surface  wave  where  c  and  V  are  the  respective  phase  speeds. 
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The   surface-wave    speed  V    is    dispersive   and   the   governing 

equation    [8]    is 

T  y 
V2   =    (g/y    +  -™   )    tanh    (V^    )  (2.l8) 

2 
where      g  =  acceleration  resulting  from  gravity  980  cm/sec 

p  =  density 

2 
T   =  surface  tension  75  dynes/cm  . 
s 

Because  of  this,  y    varies  with  Cl/2n    as  shown  in  Figure  4. 
In  this  experiment,  CI   is  of  magnitude  10  rad/sec  and 
CO  »  90  x  10   rad/sec.   Therefore,  y  ~  k  and  6d/fc2  ~  90  x  10  .  •  For 
this  reason,  the  perturbation  term  in  Eq.  2.13  appears  as  a  low- 
frequency  modulation  of  the  acoustic  signal,  and  a    becomes  a 
measure  of  the  modulation  amplitude.   The  great  inequality  in 
acoustic  and  surface  wave  speeds  means  that  the  surface  can  be 
treated  as  effectively  "frozen"  with  respect  to  the  acoustic 
pr  ocesses . 

The  solution  for  G0*  in  terms  of  £F(x,y,t)  obeys  super- 
position.  Thus,  if  the  surface  wer e  a  spectrum  of  fi,  the  solution 
for  £0   would  be  superposition  of  the  solutions  for  each  component 
of  CI   according  to  Eq.  2.13. 

In  particular,  for  a  narrow  band  of  surface  wavelengths,  a 
localized  disturbance  on  the  surface  leads  to  a  local  disturbance 
in  the  acoustic  standing  wave.   Scr imger  [5]  observed  that  the 
fluctuations  in  the  acoustic  signal  were  strongest  when  generated 
surface  waves  were  in  the  vicinity  of  the  source  or  receiver.   This 
localization  of  effect  is  assumed  to  apply  in  this  experiment. 
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III.       PROCEDURE    AND   ANALYSIS 

A.       THE   EXPERIMENT 

The    goal    of    the    experimental    investigation  was    to   achieve 
quided-mode  propagation    in   a    shallow  water    layer,    disturb   one 
surface   with   a   measurable  perturbation,    and   observe    the   resultant 
fluctuations    in    the   acoustic  pressure   amplitude. 

In   the   original    concept    of    the    experiment    design,    one    end   of   a 
shallow   tank   with  pressure-release   boundaries    was    to   be    lined  with 
sound  absorbing   mater ialand  a  fan   manifold  for    generating   wind- 
driven    "seas"   was    to   be   suspended   over    the    other    end.      By    using   a 
highly    directive   source   aimed   into    the   anechoic    area,    a    cylindri- 

cally   spreading   acoustic    traveling   wave,    which    excited   the   n      =    1 

j 

family  of  modes,  could  be  transmitted.   This  approach  would  have 
yielded  the  following  advantages: 

1.  Elimination  of  effects  resulting  from  reflections 
at  the  vertical  sides  of  the  tank. 

2.  A  simple  traveling  wave  configuration  corresponding  to 
propagation  in  the  lowest  mode  of  a  cylindr ically 
diverging  wave. 

3.  Limitation  of  the  ensonified  surface  to  a  region 
excluding  that  in  the  vicinity  of  the  fans  used  to 
generate  the  surface  fluctuations. 
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It  was  impossible  within  the  available  time  to  construct  a 
sufficiently  anechoic  tank  in  the  frequency  range  of  interest,  so 
the  method  had  to  be  abandoned.   The  alternative  was  to  place  a 
cylindr ically  symmetrical  source  in  a  corner  of  the  tank  and  ex- 
cite a  normal  mode  of  the  bounded  cavity.   This  method  raised  the 
immediate  problem  of  the  effects  of  the  highly  agitated  surface 
directly  under  the  fans.  (See  Figure  5.)   (The  surface  under  the 
fans  was  beaten  about  by  the  air  stream  generated  by  the  fans 
directed  at  the  water  surface.) 

The  theory  suggests  (as  discussed  earlier)  that  it  may  be. 
possible  to  assume  that  for  small  CF(x,y,t),  the  cause  of  the 
acoustic  amplitude  fluctuation  observed  at  the  receiver  is  the 
surface  in  the  immediate  vicinity.   Therefore,  acoustic  noise  from 
the  turbulence  under  the  fans  could  be  treated  as  a  "distant 
source"  combining  incoherently  with  the  effects  generated  by  the 
local  surface. 

Excitation  of  the  n  .=  1  family  of  modes  was  a  primary  con- 
sideration; the  complications  coming  from  model  interference  were 
to  be  avoided.  To  ensure  efficient  excitation  of  these  modes,  a 
source  active  over  a  height  greater  than  the  depth  of  the  cavity 
was  chosen.   This  insured  that,  as  the  depth  at  the  source  changed 
with  passing  waves,  source  strength  exciting  the  desired  normal 
mode  would  be  smoothly  changing  and  closely  related  to  the  local 
surface.   Since  sea  state  at  the  source  was  essentially  uncor- 
rected with  that  at  the  receiver  in  the  cavity,  and  since  the 
propagation  speed  of  the  acoustic  energy  far  exceeded  the 
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speed   of    the   surface   waves,    the   effect   was    the   equivalent    of 
simply   changing    the   amplitude   of   0    .      Equation    2.13   must    then 
be  modified   to 

$    =   [l  +  c      cos  (Q  t  +  lb    )]cos   tot    cos    k   yll+a   cos  (£2t  +  0 )  ]  (3-1) 

where   a    is    the   amplitude   fluctuation   at    the   receive   resulting   from 
the    changing    conditions    at    the   source,    and  0      and  lb    axe   assumed   in- 
coherent.      For    small    c      and    a    this    becomes 

s 


=   cos   tot    cos    k   y[l  +   \a         +    a      cos(Ot+  0')]. 


(3.2) 


In   a    similar    manner,    the    "distant   acoustic   sources"   contr-ibuted 
by    the    strong    surface-fluctuations    in    the   neighborhood   of    the    fans 
combine    incoherently   at    the   receiver    with   the  acoustic   amplitude 
fluctuations    resulting    from    the    local    surface   state.       Inclusion 
of   these    "distant    sources"   modifies    Eq .    3-2    to 


=    cos   Cut    cos    k   y[l+  \a    +    a        +    a  cos(£2t  +  t6")] 


(3.3) 


where    a  ,.   =   the   amplitude   fluctuation   at    the  receiver    resulting    from 
the    distant   sources. 

Thus,    the   sound-amplitude    fluctuations    resulting    from   the    local 
surface   state  must    be   obtained   from   the   measured  fluctuations    and 


\       2  2 

-\    a      +    ~a~      +    a 
»  r  < 


and  results    in 


-\ 


2  2 

a  -      -    a 
f  s 


(3.4) 


(3.5) 


Extraction    of    the    desired   values    of    a   and   then   e    from   the    ex- 
perimental   data    is    accomplished  by   the   following    steps: 
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1.  Measure  the  statistics  of  the  sea  surface  to  obtain 
e  as  a  function  of  surface  wave  frequency  (^2); 

2.  Measure  the  acoustic  amplitude  fluctuations  (a  )  at 

a 

the   various   positions    in    the    tank; 

3-  Determine   the    values    of    a     ,    and   a_   as    functions    of 

s        f 

surface  wave  frequency; 

4-  Compute  a   by  Eq.  3-5  and  divide  it  by  e  at  each 
frequency  of  interest. 

B.   DESCRIPTION 

The  experiment  was  conducted  in  a  shallow  water  tank  whose 
dimensions  are  shown  in  Figure  6.   The  tank  was  constructed  of 
3/4  inch  marine  plywood  set  on  2  inch  x  6  inch  supports  to  pre- 
vent sagging.   Water  depth  could  be  maintained  uniform  to  within 
jf  1  mm  throughout  the  tank. 

'A  3/8  inch  neoprene   envelope,  set  in  the  tank  and  held  in 
place  by  hydrostatic  pressure,  served  as  a  pressure  release 
boundary.   The  reflection  coefficient  of  the  neoprene  was  measured 
by  terminating  an  acoustic  waveguide  with  a  sample  of  the  material 
and  conducting  reflectivity  tests.   Values  of  around  0.95  for  the 
reflection  coefficient  were  obtained.   The  boundary  condition  in 
the  tank  was  confirmed  by  probing  the  sound  field  and  verifying 
that  the  pressure  extrapolated  to  zero,  within  experimental  accu- 
racy, as  the  probe  was  brought  near  a  boundary. 


The  neoprene  rubber  used  is  foam  containing  nitrogen  and  is 
the  type  used  in  suits  for  SCUBA  divers.   It  is  commonly  known  as 
"wetsuit"  material. 


24 


A 


117  an 


Top 


f 


237 
cm 


>- 


I 


29  cm 


Side 


Inside  Dimensions 


3/4  inch  Marine 
Plywood  


Fiberglass 
Coating 


3/8  inch  Neoprene 
-  Rubber 


Air  Gap 


Tank  Construction 
Figure  6 


25 


A  "wind-driven  sea"  was  created  by  five  equally-spaced 
squirrel-cage  blowers  in  a  manifold  at  one  end  of  the  tank.   The 
manifold  was  suspended  over  the  tank  to  prevent  coupling  of 
mechanical  vibration  to  the  tank  (Figure  7).   Each  fan  had  a  rated 
capacity  of  100  ft  /min.   Various  seas  could  be  generated  by  varying 
the  angle  of  attack  of  the  fan  output,  the  height  of  the  manifold, 
orifice  size,  and  the  number  of  fans  (Figure  8).   Greater  sea 
heights  were  obtained  by  placing  a  16  inch  wide  cover  over  the 
water  adjacent  to  the  fan  manifold  in  order  to  hold  the  wind  down 
on  the  surface  (Figure  5) • 

The  wave  heights  were  measured  by  an  immersion  probe  whose 
impedance  varied  inversely  with  its  submerged  length.   This  wave 
probe  consisted  of  two  submerged  electrodes  with  a  10  kHz  cur- 
rent flowing  between  them.   The  passing  surface  waves  modulated 
the  current  amplitude  by  changing  the  amount  of  submerged  elec- 
trode.  This  produced  a  current  whose  envelope  amplitude  was  a 
measure  of  the  passing  wave  height.   After  amplification,  the 
information  was  retrieved  by  filtering  and  detecting  the  amplitude 
envelope  of  the  current. 

Wave  probe  construction  is  shown  in  Figure  9.  The  positive 
electrode  was  made  from  20-mil  nikomel  wire  soldered  to  1/16  inch 
brazing  rod.  The  gound  return  was  1/16  inch  brazing  rod  coated 
with  neoprene  paint  over  its  vertical  length.  The  last  4  cm  were 
kept  bare  and  bent  at  a  right  angle  to  the  vertical.  The  elec- 
trodes were  soldered  to  standard  banana  plugs  and  held  in  place 
by  a  connector  (Pomona  Electronics  Connector  model  #l68l) . 
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Tank  with  Fan  Manifold  in  Place 
Figure  7 
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The  Wave  Probe 
Figure  9 
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The    detection   and  filtering    network    (see   Figure    10)    con- 
sisted  of   a   switching    diode   for    half -wave   rectification,    two 
sucessive    low-pass    filters    with    cutoff   frequencies    of    159.5   Hz, 
and  a    blocking   capacitor.      The   25  M<F  capacitor    with    the    50    k  Cl 
input    impedance    of    the   following   amplifier    formed  a    high-pass 
filter    with  0.126    Hz    cutoff.      The   network  was    tested  with   an 
unmodulated   10   kHz    tone  which    ideally   would  generate   zero    DC  and 
AC   outputs.      The   AC   signal   which    leaked   through   the    filter    amounted 
to  no   more    than   0.05  percent    of   the    input   amplitude,    and   the   DC 
appearing   at    the   output   was    equally    insignificant. 

The   relation   between    submerged    length    of    the  wave   probe  and 
current    flow    was    non-linear,    the   effect   being   more    or    less   severe 
depending    on    the  positioning    of    the  probe   with   respect    to   the 
cavity    boundaries.       (If   the   electrodes   were  placed   too   near    the 
bottom,    the   non-linearity   was    aggravated.)       By   bending    the    ground 
return   and  adjusting   the  probe   current    level,    an   adequately    linear- 
operating   range   was   established   for    the   wave   heights    expected. 

Due   to   the   bottom   effects    mentioned,    a    static   probe   cali- 
bration  method  was   used:    With    the  probe    in  place,    the   water    level 
in    the    tank  was    varied  and   the    output    voltage   measured   on    the 
probe   side    of   the    blocking   capacitor.       (It    was    verified  by   use    of 
an  AN/URM-26   AM  modulation    test    set    that    there  was    no  measurable 
loss    in    amplitude  across    the   blocking    capacitor    in    the    frequency 
range   of    interest.) 
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LIST    OF   EQUIPMENT 
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Amp  .  2 
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0.1  Octave  Filter 

XY  Plotter 
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WAVETEK  Multipurpose  VCG  Model  116 

HEWLETT  PACKARD  2470A 

HEWLETT  PACKARD  46 7A 

TEKTRONICS  Model  1121 

HEWLETT  PACKARD  46 OA 

PRINCETON  APPLIED  RESEARCH  Model  100A 

PRECISION  INSTRUMENT  6200 

AMPEX  PR-10 

GENERAL    RADIO  Model    1564A 

MOSELEY   7035B 

HEWLETT   PACKARD  120B 

ALLISON  Model    2C 

HEWLETT   PACKARD   427A 

HEWLETT   PACKARD   3400A 
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The  system  for  processing  wave  information  is  shown  in 
Figure  11.   The  wave  information,  contained  in  the  envelope  of  the 
wave  probe  signal,  was  first  amplified;  this  information  was  then 
recorded  on  an  FM  tape  recorder  and  simultaneously  sent  to  a  signal 
correlator.   The  output  of  the  correlator  was  graphically  recorded 
on  the  XY  plotter. 

The  sound  source  used  was  a  BaTiO   cylinder,  4-5  cm  in 

diameter  and  10.5  cm  tall.   The  receiver  was  a  1/8  inch  BaTiO 

3 

probe.   The  sound  processing  system  is  shown  in  Figure  12.   The 
received  signal  was  pre-amplif ied  in  a  low-noise  40  db  amplifier, 
bandpass  filtered  at  carrier  +  2.5  kHz  and  then  amplifier  further. 
The  amplitude  modulation  of  the  acoustic  signal  was  analyzed  in 
like  manner  with  the  same  detector  network  used  in  the  wave  probe 
system.   This  information  was  then  recorded  in  the  same  manner  as 
the  wave  height  data. 

The  wave  and  acoustic  amplitude  fluctuation  data  were  analyzed 
for  frequency  content.   The  system  used  is  shown  in  Figure  13.   The 
original  tapes,  recorded  in  the  systems  previously  described,  were 
replayed  100  times  faster  and  recorded  on  tape  loops  with  a 
standard  audio  recorder.   The  tape  loops  were  then  replayed  through 
a  0.1  octave  filter  as  shown. 

C.   PROCEDURE 

The  experiment  was  conducted  in  water  of  8.8  cm  depth.   A 

strong  normal  mode  was  found  at  l4.1  kHz,  a  frequency  well  below 

cutoff  of  the  n   =  2  family  of  modes.   The  sound  field  was  probed 
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in  order  to  identify  the  modes  nodel  pressures  other  than  zero 
were  found  in  the  tank,  indicating  excitation  of  more  than  one 
mode  or  more  likely  imperfections  in  the  boundaries.  The  cri- 
terion for  a  strong  mode  was  taken  to  be  the  average  standing  wave 
ratio,  which  in  this  case  was  approximately  60. 

The  mode  was  identified  by  measuring  the  lengths  of  the 

standing  waves  and  estimating  the  wave  numbers.   The  mode  with 

n   =  28,  n   =  95  n   =  1  appeared  to  be  excited.   This  mode  has  a 
1         2        3 

calculated  frequency  of  13. 6  kHz  (4*25  percent  less  than  the 
observed  frequency).   The  vertical  sound  profile  ver if ied #that 
the  n   =  1  family  of  modes  was  being  excited  (Figure  14)  as 
assumed  in  Section  II. 

Once  the  desired  mode  was  excited  and  identified,  three  re- 
ceiver stations  were  selected  at  strong  antinodes  of  the  hori- 
zontal wave  pattern  (Figure  15).   The  procedure  now  described  was 
followed  at  each  receiver  position. 

The  acoustic  probe  was  set  at  a  depth  of  l\.i\    cm  in  still  water, 
and  the  sound  carrier  amplitude  measured.   The  surface  was  then 
disturbed  with  a  preselected  wind  condition.   The  following  ob- 
servations were  then  made  of  the  envelope  of  the  modulated  sound 
carrier : 

1.   A  500  second  record  was  made  of  the  mean  square 

amplitude  of  the  modulation  from  the  output  of  the 
correlator  on  the  XY  plotter. 
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2.  A  20  minute  FM  tape  recording  was  made  of  the 
acoustic  modulation. 

3.  The  autocorrelation  function  was  computed  at  two 
different  times  by  the  correlator  and  recorded 
by  the  XY  plotter . 

The  wind  generator  was  then  stopped  and  the  surface  allowed 
to  calm.   The  carrier  amplitude  was  again  measured  to  verify  that 
the  amplitude  of  the  standing  wave  had  not  changed.   The  acoustic 
data  for  six  different  sea  states  at  each  receiver  position  were 
obtained  in  this  fashion. 

The  area  directly  under  the  fans  of  the  wind  generator  was 
intensely  agitated  relative  to  the  rest  of  the  sea  surface.   In 
order  to  isolate  this  contribution  to  the  modulation  of  the  signal, 
the  tank  was  modified  as  shown  in  Figure  16 .   The  reverse  flap 
lifted  the  wind  off  the  surface,  and  the  styrofoam  prevented  the 
propagation  of  surface  waves.   Since  the  styrofoam  was  pressure- 
release,  its  presence  had* a  negligible  effect  on  the  standing  wave 
pattern . 

With  the  probe  in  position  at  a  designated  receiver  station 
and  the  wind  machine  operating  at  a  specified  setting,  the  ampli- 
tude of  the  acoustic  modulation  was  measured  through  a  0.1  octave 
filter  and  average  modulation  amplitude  as  a  function  of  frequency 
was  recorded.   These  data  were  recorded  at  positions  "B"  and  MC" 
only;  position  "A"  was  obstructed  by  the  styrofoam  buffer. 


40 


kl 


Once  all  sound  data  had  been  collected,  the  statistics  of  the 
sea  surface  were  measured.   A  calibration  of  the  wave  probe  was 
always  conducted  before  a  data  run.   With  the  probe  in  place,  the 
tank  was  filled  with  water  to  a  depth  of  10  era.   The  water  was 
then  slowly  siphoned  out  and  the  DC  output  of  the  wave  probe  was 
recorded  as  a  function  of  depth.   A  calibration  factor  was  de- 
termined from  this  curve.   Probe  performance  did  not  vary  with 
position  in  the  tank  except  very  near  the  walls. 

A  static  calibration  of  the  wave  probe  is  considered  valid; 
however,  a  dynamic  check  was  performed  by  simultaneously  measuring 
the  height  of  a  passing  wave  visually  and  comparing  it  with  the  . 
output  of  the  probe.   The  passing  waves  were  measured  on  a  metric 
scale,  printed  on  heavy  photographic  paper  and  taped  to  the  side 
of  the  probe  holder.   (The  thin  paper  ruler  caused  negligible  dis- 
turbance in  the  passing  surface  waves.)   The  output  of  the  probe 
was  observed  on  an  oscilloscope. 

Once  the  probe  "Was  calibrated,  it  was  placed  at  each 
receiver  position  and  measurements  were  made  of  the  amplitudes  of 
the  waves  generated  on  the  surface. 

The  statistics  of  the  surface  at  the  position  of  the  sound 
source  were  also  required.   These  data  could  not  be  taken  with 
the  wave  probe,  so  the  behavior  of  the  various  frequency  components 
of  the  sea  were  measured  as  they  propagated  down  the  tank  and  a 
projection  of  the  behavior  of  each  water  surface  at  the  position 
of  the  acoustic  source  was  made.   This  estimate  was  checked 
visually  and  appeared  to  be  valid  within  experimental  error. 
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Analog  frequency  analysis  of  the  FM  tape  recording  was  made 
by  a  tape  loop  method.   The  20  minute  tape  records  were  played 
back  100  times  faster  and  recorded  on  7  second  tape  loops  with  and 
audio-magnetic  tape  recorder. 

These  loops  were  then  played  through  a  0.1  octave  filter  with 
slow  integration  time.   Average  amplitude  as  a  function  of  fre- 
quency was  recorded.   The  tape  loop  output  was  also  played  through 
a  true  RMS  voltmeter  to  measure  overall  RMS  amplitude  for  com- 
parison with  the  mean  square  amplitude  prediction  of  the  correlator 

D.   ANALYSIS 

The  RMS  wave  height  was  computed  by  having  the  correlator  set 
to  compute  the  autocorrelation  function.   In  this  mode  the  cor- 
relator output  is 

T 

c(T)  =  i  f   V.n2(t)dt  (3.6) 

J   o 

which  equals  the  mean  square  value  of  the  input  when  T  approaches 
infinity.   The  correlator  integration  time  is  only  20  seconds,  so 
a  good  estimate  of  the  mean  square  value  was  obtained  by  taking 
the  average  of  a  500  second  graphical  record  of  the  correlator 
output  (Figure  17).   The  various  equipment  amplification  factors 
were  divided  out  and  the  square  root  taken. 

The  RMS  value  of  the  envelope  amplitude  of  the  wave  probe 
signal  was  computed  in  the  same  manner. 

The  average  surface  wave  and  acoustic  amplitudes  as  functions 
of  frequency  were  obtained  as  previously  discussed.   As  the  data 
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were  recorded,  a  100  Hz  tone  was  recorded  on  another  channel  of 
the  tape.   The  system  amplification  was  checked  by  monitoring  this 
test  signal  at  various  stages  in  the  frequency  analysis.   This 
check  also  giwe   assurance  that  the  data  were  not  being  distorted 
during  the  processing. 

In  order  that  a   and  C  be  computed,  the  gains  of  the  equipment 
through  which  the  data  was  sent  had  to  be  divided  out  of  the 
observations . 

Refer  to  Figures  11  and  12  and  let 

GAINA„,„  =  the  amplifier  gain 
AMP 

GAIN       =  the  FM  tape  recorder 
lAPb  -L 

GAIN^A„_  _  =  the  audio  tape  recorder. 
TAPE  2 

The  average  amplitude  measured  through  the  0.1  octave  filter 

is  converted  to  the  average  value  of  e  for  each  frequency  interval 

as  follows: 

WAVE  PROBE  OUTPUT  =  FILTER  0UTPUT 


GAINAMP  •  GAINTAPE  1  •  GAINTAPE  2 


The  wave  probe  output  can  be  expressed  in  centimeters  of 
wave  height  by  dividing  by  the  DC  calibration  factor  determined 
before  each  run: 

WAVE  HEIGHT  =  -S^^-P^^I 

AV  CALIBRATION  FACTOR  * 

The  normalized  fluctuation  e  then  becomes 

\  /   l \  /: 


e  = 


[NAMP-    GAINTAPE    T    GAINTAPE    z){ CAL   ™CT0RJpPTH 

(3-7) 
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FILTER    OUTPUT 


^"W    CAL    FACT0R     '     DEPTH      VGAINTAPE    1    GAINTAPE   2) 

The   sound  amplitude   fluctuation    a   at    each   frequency    of 

interest    is 

_  ENVELOPE  AMPLITUDE  AT  THAT  FREQUENCY 
a  CARRIER  AMPLITUDE 

Since  the  sound  information  was  processed  with  the  same  system 

as  the  surface  data,  a   was  computed  as  follows: 

_  AVERAGE  AMPLITUDE         1 g 

'  '    <GAINAMP><CARRAMP>    (GAINTAPE  1»(GAINTAPE2'  "  '    \ 

Notice  that  no  calibration  factor  is  required  in  the  calculation 

of  a  . 

In  the  ratio  a/e  ,  the  factor  -—•■■; Trrrr, vanishes; 

GAINTAPE  1  GAINTAPE  2 

therefore,  all  data  were  treated  in  arbitrary  units  with  this 

factor  included. 

The  carrier  amplitude  by  which  the  envelope  amplitude  was 
divided  was  an  average  of 'the  amplitude  measured  before  and  after 
each  run.   If,  however,  the  difference  in  these  amplitudes  was 
greater  than  2  percent,  the  run  was  re-executed. 

The  presence  of  high  frequency  sound  fluctuations,  unpre- 
dicted  by  the  theory,  and  the  persistence  of  these  fluctuations 
throughout  the  tank  strongly  suggested  the  fans  as  a  source.   The 
acoustic  amplitude  fluctuation  due  to  the  agitation  under  the  fans 
was  obtained  for  positions  "B"  and  "C"  and  described  above.   For 
any  given  sea  state,  these  data  were  very  similar  at  the  two 
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positions.      Therefore,    an   average    (fan    noise)    level    for    each   sea 
state   was    assumed    (Figure    18)    for    all    the    receiver    stations. 

The   other    factor    to   be  accounted   for    was    source  modulation. 
The   sea    state   at    the   source  was   predicted  as   previously    described.. 
The    effect    of    the   waves    on    the   acoustic    output    of    the   source  was 
determined  by    raising    the    source   a    given   distance    out    of    the  water 
and  measuring    the   reduction    in   sound  amplitude   at  position    "B". 
With   a    10   percent   reduction    in    the    radiating    surface,    the    sound 
amplitude  at   position    MB"    was    reduced  by   6   percent.      This    relation 
was    found   to   be    linear    within    the   range    of    interest.      Therefore, 
the   amplitude   modulation    of    the    source  was   assumed   to   be 

as   =  0.6   c,  (3.9) 

(This  result  was  checked  by  shading  the  upper  part  of  the 
transducer  the  same  amount  with  3/8  inch  neoprene  foam  rubber  and 
comparing  the  results,  which  agreed  within  experimental  error.) 
The  amplitude  fluctuation  due  to  the  local  surface  effect  was  then 
computed  with  Eq .  3*5: 


''ca   "  cs   "  af  (3-5) 

As  a  check  on  the  validity  of  the  data,  several  quantities 
were  measured  by  two  methods  and  the  results  compared.   The  tape 
loops,  for  instance,  were  played  through  the  true  RMS  voltmeter  in 
order  to  measure  the  RMS  amplitude  fluctuation  and  wave  height. 
The  results,  adjusted  for  system  gain,  agreed  to  within  10  percent 
with  the  RMS  values  computed  through  the  correlator. 
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The  frequency  distribution  of  the  modulation  data  was  also 
determined  from  the  autocorrelation  function  plotted  during  each 
data  run.   The  cosine  Fourier  transform  of  the  autocorrelation 
function  c(T)  equals  the  power  spectral  density.   The  transform 
of  c(T)  was  computed  by  numerical  integration  of  100  points  taken 
from  the  plots  described  above.   The  resulting  spectra  confirmed 
the  distributions  obtained  from  the  tape  loops. 

E.   ERROR  ANALYSIS 

The  major  cause  of  error-  in  this  experiment  was  the  non- 
stationarity  of  the  processes.  The  variation  in  mean  square 
amplitude  (see  Figure  17)  gives  a  qualitative  measure  of  the 
degree  of  non-stat ionar ity  involved. 

Measuring  equipment  with  sufficiently  long  integration  times 
to  overcome  this  problem  was  not  available.   Therefore,  long 
"effective"  integration  times  were  achieved  by  visual  integration 
of  the  meter  readings  or  .graphical  averaging  as  previously  dis- 
cussed.  Study  of  Figure  17  suggests  that  an  integration  time  of 
500  seconds  is  sufficient  for  our  purposes. 

In  the  frequency  analysis,  15  minutes  of  FM  tape  record  was 

compressed  onto  a  7  second  tape  loop  (see  page  33)-   Visually 

averaging  the  meter  reading  on  the  0.1  octave  filter,  provided  an 

integration  time  in  excess  of  500  seconds.   The  estimated  error  in 

the  value  assigned  to  acoustic  amplitude  fluctuation  ( a    )  or  sur- 

a 

face  height  fluctuation  (e)  is  _+  10  percent.   (This  estimate  was 
determined  by  comparing  the  values  estimated  by  two  independent 
observers . ) 
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The  tape  loop  method  was  not  used  in  obtaining  the  acoustic 
amplitude  fluctuation  due  to  the  fans;  in  this  case  the  filter  was 
placed  directly  on  the  output  of  the  filtering  and  detecting 
network  (page  40)  .   Averaging  two  sets  of  observations  together 
at  each  wind  setting,  minimized  the  error  in  the  fan  contribution 
(a_p)    to  an  estimated  15  percent,  so  that  the  more  elaborate  tape 
loop  analysis  was  not  deemed  necessary. 

Recall  that  the  acoustic  amplitude  fluctuation  term  a 
resulted  from  the  source  modulation,  which  was  extracted  from 
the  e  curves.   The  estimated  error  in  £  is  _+  10  percent.  'Adding 
a  _+  5  percent  error  for  the  graphical  extrapolation  yielding  the 
surface  at  the  source  position,  results  in  a  total  error  in  the 
source  modulation  term,  a    ,  of  +  15  percent. 

The  RMS  amplitude  values  computed  from  the  graphical  output 

of  the  correlator  are  considered  reliable.   In  order  to  test  the 

validity  of  the  average  amplitudes  obtained  from  the  tape  loops, 

RMS  amplitudes  computed  in  this  manner  were  compared  to  those 

from  the  correlator.   The  values  were  consistently  within 

10  percent . 

The   error    in    a    ,    a    ,    a^.   and  e    were   assumed   statistically    in- 
a'      s         f 

dependent   at    each   receiver    position.      Therefore   the   error    in    the 
value    of   a/e    at    each   frequency   was    computed  using    standard  RMS 
error    analysis . 
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IV.   RESULTS  AND  CONCLUSIONS 

The  experimental  results  consist  of  eighteen  sets  of  data:    , 
measurements  were  taken  at  three  different  receiver  positions  for 
six  different  winds.   The  six  winds  generated  can  be  grouped  in 
pairs:  (1,2),  (3,4)  and  (5,6).   Each  pair  of  winds  was  generated 
with  the  same  settings  of  manifold  height,  angle  of  attack,  and 
orifice  size  of  the  wind  machine.   The  even-numbered  winds  were 
generated  with  five  fans;  the  odd-numbered  ones,  with  three.   The 
receiver  positions  were  previously  shown  in  Figure  15- 

The  frequency  distribution  of  the  measured  acoustic-amplitude 

fluctuation  (a    ),  the  fluctuation  induced  by  the  fans  (a.-),  and  the 
a  x 

amplitude  fluctuation  due  to  source  modulation  (a    )  are  plotted 
for  each  wind  setting  at  a  given  receiver  position  (see  Appendix  P>)  . 
These  graphs  reveal  certain  general  features: 

1.  The  frequency  spectrum  of  the  measured  acoustic 

fluctuation  (a    )    closely  resembles  that  of  the  local 
v  a' 

surface  fluctuation  (c).   In  all  cases  each  curve  computed 
for  a   closely  resembles  the  curve  for  the  equivalent  ar   • 

2.  At  each  receiver  position,  the  frequency  spectra  for 
both  "B"  and  "C"  are  similar  for  each  wind  pair.   The 
stronger  winds  produce  greater  amplitudes  as  would 

be  expected. 

3.  The  curves  of  a      are  similar  at  all  receiver  positions 

a  f 

for  any  given  wind  setting.   The  difference  in  values 

of  a      for  positions  "B"  and  "C"  is  very  small. 
a      ^  * 
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4.  The  higher  frequency  components  of  e  appear  to  die  out 

as  the  surface  waves  progress  down  the  tank.   (See 

Figure  22.)   The  corresponding  values  of  a      fall  off 

a 

much  more  slowly. 

5.  In  computing  the  values  of  the  acoustic  amplitude 
fluctuation  resulting  from  the  local  surface 


\(a)    -    (a    )    -  (a    f    ,  (3-5) 


'f 

we  see  that  a   becomes  significant  at  the  higher 
frequencies  while  a      is  significant  at  the  lower 
ft equencies . 
6.   There  appears  to  be  evidence  of  generation  in  a      of 
second  harmonics  of  frequencies  in  the  region  corre- 
sponding to  resonance. 
The  computed  values  of  a/e  (using  Eq.  2.13)  are  plotted  as 
functions  of  surface-wave  frequency  for  each  wind  condition  at 
each  receiver  position.   (Appendix  C.)   The  results  are  summarized 
in  Figure  19  where  the  curves  are  collected  on  single  graphs 
according  to  wind  condition. 

In  the  region  below  the  predicted  "resonance"  the  measured 
values  of  a/e  compared  favorably  with  the  theoretical  prediction. 
Very  little  evidence  for  a  resonance  effect  exists  in  the  results 
of  any  run.  Above  the  predicted  resonance  the  observed  values  of 
a/e  are  consistently  greater  than  the  theoretical  values.  The 
absence  of  any  peaking  of  a/e  in  the  vicinity  of  resonance  cannot 
be  explained. 
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Several    important    factors    could  be   contributing    to    the   high- 
frequency    behavior    mentioned.       Recall    that    the    term    o       is    ex- 
tremely   important    in    determining    the   acoustic   amplitude    fluctu- 
ation   a    in    this   region.       It    is    suspected   that   perhaps    the   entire 

observed  acoustic   amplitude   fluctuation    a      above   6    Hz    can   be 

^  a 

attributed  to  the  surface  agitation  under  the  fans. 

In  the  design  of  the  fan  isolation  experiment,  it  was  assumed 
that  the  surface  agitation  directly  under  the  fans  was  responsible 
for  the  additional  amplitude  fluctuation  at  the  receiver.   The 
styrofoam  buffer  was  therefore  placed  under  the  orifice  flap. 
This,  however,  suppressed  large  surface  waves  generated  in  this 
region.   It  is  conceivable  that  these  waves  contribute  to  the 
"distant  sources",  but  this  possibility  was  not  investigated  at 
the  time.   Recall  also  that  analyses  of  the  fan-generated  noise 
were  not  conducted  with  the  more  accurate  tape- loop  method 
because  of  time  and  equipment  limitations.   Therefore,  the  a~ 
term  may  be  larger  than  recorded.   This  hypothesis  would  have 
been  investigated  further  if  time  had  been  available. 

Since  accuracy  in  a,   becomes  critical  when  computing  the 
ratio  c/e  for  the  higher  frequencies,  increasingly  anomolous 
behavior  in  a      in  this  frequency  range  at  distances  further 
removed  from  the  fans  cannot  be  considered  as  unequivocably  con- 
tradicting the  theory. 

The  discrepancy  between  experiment  and  theory  for  values  of 
c/e  at  the  higher  frequencies  could  also  be  attributed  to 
violations  of  the  assumptions  in  the  theory.   For  instance, 
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Figure  20  shows  that  the  surface  fluctuations  were  functions  of 
both  x  and  y,  rather  than  just  of  x  as  was  assumed.   Investigation 
of  a  more  realistic  expression  for  F(x,.y,t)  would  require  extensive 
auto-  and  cross-correlation  on  measurements  of  the  generated  sur- 
face fluctuation. 

It  is  also  assumed  in  the  theory  that  the  properties  of  the 
surface  fluctuations  are  uniform  in  the  vicinity  of  the  receiver. 
Figure  21  shows  the  distribution  of  RMS  wave  heights  for  a  typical 
wind  over  the  area  of  the  surface  in  the  vicinity  of  the  receiver 
stations.   This  figure  shows  that  the  waves  are  dying  out  as  they 
propagate  down  the  tank.   Figure  22  shows  that  the  frequency 
spectrum  for  the  surface  waves  above  the  frequency  corresponding 
to  the  "resonance"  is  also  changing.   The  effect  of  the  inhomo- 
geneity  of  the  surface  could  not  be  determined  but  is  probably 
most  severe  at  position  "A". 

A  basic  assumption  in  the  theory  is  that  the  mechanism  for  the 
acoustic  fluctuation  is  the  local  surface  state.   In  the  investi- 
gation of  the  validity  of  this  assumption,  the  following  points 
merit  study: 

1.   With  the  normal  mode  excited  and  the  water  surface 

calm,  the  surface  in  the  vicinity  of  the  source  was 

2 
disturbed  by  quadradigital  agitation  .   An  immediate 

amplitude  fluctuation  was  observed  at  the  distant 


2 
Trickling  four  fingers  in  the  water,  making  small  capillary 


waves  . 
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receiver  (station  "B").   The  surface  above  the  receiver 
was  disturbed  in  the  same  manner  and  the  acoustic  signal 
was  again  immediately  perturbed.   However,  when  the  water 
surface  was  disturbed  at  a  point  removed  from  both  source 
and  receiver,  an  effect  on  the  received  acoustic  signal 
was  observed  only  after  the  surface  disturbance  had 
propagated  into  the  vicinity  of  the  source,  receiver, 
and  region  in  between. 
2.   The  treatment  of  the  extreme  agitation  under  the  fans 
as  a  distant  source  (page  lh)    contributing  to  observed 
acoustic-amplitude  fluctuation  appears  valid  from  the 
results  of  the  fan  isolation  experiment,  but  this  would 
appear  to  be  inconsistent  with  the  result  just  discussed. 
Recall,  however,  that  the  violent  agitation  of  the  sur- 
face under  the  fans  is  a  clear  violation  of  the  assumption 
that  CF(x,y,t)  be  small. 
The  first  point  above  appears  to  validate  the  assumption  that 
the  driving  mechanism  for  the  acoustic  fluctuation  at  the  receiver 
is  the  local  surface.   The  fan  isolation  experiment- (point  two) 
showed  that  severe  surface  disturbances  at  a  remote  location  could 
cause  amplitude  fluctuations  under  a  smooth  surface.   It  is  plau- 
sible that  the  distinction  between  these  two  cases  results  from 
the  difference  in  magnitude  of  GF(x,y,t)  in  the  two  regions. 

The  stringent  assumptions  of  the  simple  theory  were  all  clearly 
violated  to  some  degree  during  the  conduct  of  the  experiment.   It 
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is  concluded,  however,  that  the  measured  values  of  a/e  below  the 
frequency  of  predicted  "resonance"  are  sufficiently  close  to  the 
theoretical  values  to  claim  consistency  between  experiment  and 
theory  in  this  range. 


65 


APPENDIX  A 
Computer    Program   for    a/e    Calculation 

1.  IIPLICIT    REAL*8(A,    K,    P,) 

2.  REAL*8(DCOS,  DCOSH ,  DABS,  DSQRT) 

3.  DO  100  1=1,  500 

4.  GAM  =  (I-l)*.01 

5.  Kl=2.*3. 1415926535/16.  4 

6.  K3=3- 14159265535/8.8 

7.  PI^3- 1415926535 

8.  AR1=2.*GAM*K1 

9.  AR2=K3**2 

10.  AR3=GAM**2 

11.  AR4=Pl/2. 

12.  AR5=(1.-((AR3+AR1)/AR2)) 

13.  AR6=(1.-((AR3-AR1)/AR2)) 

14.  IF(AR5.LT.O. )GO  TO  50 

15.  AR7=DCOS(AR4*DSQRT(AR5)) 

16.  GO  TO  60 

17.  50  AR7=DCOSH(AR4*DSQRT(DABS(AR5) ) ) 

18.  60  IF(AR6.LT.O.)  go  to  70 

19.  ar8=dcos(ar4*dsqrt(ar6) ) 

20.  70  ar8=dcosh(ar4*dsqrt(dabs(ar6) ) ) 

21.  80  ar9=dabs(1./ar7-1./ar8) 

22.  ABE=PI*AR9/4. 

23.  PRINT,  GAM,  ABE 

24.  100  CONTINUE 

25.  STOP 

26 .  END 
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APPENDIX  B 
Observed  Data 
The  eighteen  sets  of  data  described  on  page  51  are  plotted  as 
functions  of  surface  wave  frequency  (Cl/2rt)  . 

a   =  the  measured  acoustic  amplitude  fluctuation 
e  =  the  normalized  surface  wave  fluctuation  in  the 
vicinity  of  the  receiver 
a~   -    the  acoustic  amplitude  fluctuation  generated  by  the 

severe  agitation  under  the  fans 
a      -    the  acoustic  amplitude  fluctuation  due  to  source 
modulation 
Representative  error  flags  for  these  data  are  presented  in 
Figure  Bl . 
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APPENDIX  C 
Measured  Values  of  a/e 

The  measured  values  of  a/e  at  each  receiver  position  for 
each  wind  condition  are  plotted  as  functions  of  surface  wave 
frequency  (Q/2n).   Standard  RMS  error  flags  are  presented. 
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